Full counting statistics of chiral Luttinger liquids with impurities 
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We study the statistics of charge transfer through an impurity in a chiral Luttinger liquid (realized 
experimentally as a quantum point contact in a fractional quantum Hall edge state device). Taking 
advantage of the integrability we present a procedure for obtaining the cumulant generating function 
of the probability distribution to transfer a fixed amount of charge through the constriction. Using 
this approach we analyze in detail the behaviour of the third cumulant C3 as a function of applied 
voltage, temperature and barrier height. We predict that C3 can be used to measure the fractional 
charge at temperatures, which are several orders of magnitude higher than those needed to extract 
the fractional charge from the measurement of the second cumulant. Moreover, we identify the 
component of C3, which carries the information about the fractional charge. 
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It was realized only recently that the probability distri- 
bution function P{Q) of transmitting a charge Q through 
a given system per unit time (also referred to as full 
counting statistics, FCS) is the ultimate transport char- 
acteristics of a given device 0. Along with the current 
voltage relation it contains all information about the fluc- 
tuations of current and transmitted charge. P{Q) of a 
clean wire is purely Gaussian, i. e. that only the first two 
cumulants (irreducible moments) are different from zero. 
On the contrary, as soon as a scatterer is inserted, the 
transport is described by a binomial distribution func- 
tion, all cumulants of which arc known to be nonzero 
P, y| • The issue of the third cumulant C3 appears to be 
especially interesting in systems coupled to the electro- 
static environment 0, Q . Crucial parameters concerning 
the influence of the latter can be investigated by the mea- 
surement of C3. 

Another question the FCS is supposed to answer is that 
of the charge of current carrying excitations. The first at- 
tempt goes back to Schottky, who already 1918 has stated 
that the proportionality factor between the current au- 
tocorrelation spectra at zero frequency and the charge 
current is the elementary charge provided the transmis- 
sion is weak j^. Using the extrapolation of his ideas 

0, to transport phenomena in fractional quantum 
Hall (FQH) bars it has been established, that the charge 
of current carrying excitations is indeed a fraction of the 
elementary charge 0,^3- Such kind of analysis is, how- 
ever, quite difficult as one has to subtract from the cor- 
relation data the Johnson-Nyquist (JN) contribution re- 
sponsible for the equilibrium thermal fiuctuations. On 
the other hand, it has been suggested in ^ by calcu- 
lations in the tunneling regime, that the third cumulant 
under similar conditions may offer much more reliable 
means of measuring the charge because it does not have 
the JN constituent. Thus far no experimental attempt to 
extract the fractional charge from higher cumulants has 
been undertaken however. 

The fractional charge measured in |l[l3 undoubtedly 
comes about as a result of electronic correlations. Until 



now the analysis of fully interacting systems has been 
restricted only to some special cases 

[S El Q III IS 

IT^ ITsI l . One way to make progress in understanding the 
FCS of fully interacting systems is to take advantage of 
integrability. Especially in ID there are quite a number 
of exactly solvable models, one notable example being 
the chiral Luttinger liquid (CLL) with an impurity, which 
describes in particular tunnelling between FQH edges 0, 
I20I I . In this Letter we attempt to use the integrability of 
this model in order to access its FCS. 

We start with the following CLL Hamiltonian in the 
bulk, which can be written in terms of right/left {R,L) 
moving charge current densities. 
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(1) 



where we fix h — 1 and where g is the Luttinger liquid 
interaction parameter In the case of a FQH device 
g = V when 1/v is an odd integer [T^ l2?il |. The origi- 
nal Fermions can be rewritten as exponents of the phase 
fields jL,R = T9i:</'L,i?/27r, so that the (local at x = 0) 
backscattering term (or interedge tunnelling) is given by 



Hbs = \b cos [M^)-<pRm 



(2) 



It has been shown in |22| , that + Hbs can be mapped 
(after a folding and a transformation into an even-odd 
basis) onto the boundary sine-Gordon model, which is 
integrable by means of Bethe ansatz (BA) technique 
[23| . At 1/v odd integer - to which we restrict now - 
the corresponding excitations are (anti) kinks (denoted 
by subscripts ±), propagating in two different directions, 
and i ~ 1 . . ~ 2 breathers. All of them are mass- 
less and have the dispersion relations e± = e^/2 and 
Cj = sm^K i / 2{1 / V — 1)] e^, where 9 is referred to as ra- 
pidity and parameterizes the particle momenta. While 
the breathers are neutral, the (anti)kinks carry charge 
(in units of the electron charge) g = =pl. Furthermore, 
similarly to the noninteracting systems one can introduce 
the level densities (per unit length and rapidity) ni{0) 
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and filling fractions fi{0), i = ±,l,...l/i/ — 2, which 
can be combined to give the density of occupied states 
Pi{9) = ni{9)fi(9). The filling fractions can, in turn, be 
rewritten in terms of quasienergies according to the 
prescription fi = Pi/rii = 1/(1 + e^'^'"'''^/-^), where 
is the chemical potential and T is the temperature. The 
equilibrium values of quasienergies as well as those of 
density of occupied states arc found by means of thermo- 
dynamic Bethe ansatz (TEA) as derived in 0| • The 
charge current / can then be shown to be carried by sin- 
gle (anti)kinks scattering processes across the boundary 
(which account for charge transfer accross the impurity 
in the initial problem), breathers only role being in es- 
tablishing the equilibrium values of P± = {P±), 



I = i{v,Tb) = I dOTio) [p+{e) - p_(0)] , 



(3) 



which arc found using TBA [H |2l III. V denotes 
the voltage applied across the constriction. The quan- 
tity T{9) is the probability for a kink or anti-kink with 
rapidity 9 to scatter off the boundary as an anti-kink or 
kink (so the current is conserved) and is given in .22j • It 
strongly depends on Tb , which is the energy scale gener- 
ated by the backscattering strength A b ■ The first cumu- 
lant of transmitted charge is then Ci = rl/h , r being 
the waiting time. The zero frequency current autocorre- 
lation spectra (sometimes also called noise) was obtained 
in p5l |. the corresponding cumulant of charge given by 
Eq. (27) of 12^, supplemented with the prcfactor r/h. 

For general systems, there are two ways to obtain the 
full P{Q)- One possibility is the calculation of the cu- 
mulant generating function x(A) = J2q e^^^P{Q), which 
can be calculated either directly via some functional in- 
tegral [23, or it can be reduced to a calculation of some 
nonequilibrium Green's function |l3l IT^ . The other way 
is a direct evaluation of time-dependent current corre- 
lation functions flj and corresponding cumulants with a 
subsequent reconstruction of P{Q), see e. g. In 
both cases one has to solve for nonequilibrium quanti- 
ties. The caveat of the first approach is the fact that 
one has to work with some kind of nonequilibrium tech- 
nique directly, e. g. Keldysh diagrammatics. Even for 
integrable systems this cannot yet be done as it requires 
the knowledge of explicitly time dependent Green's func- 
tions, which are very difficult to extract from the BA. In 
the second situation the difficulty is the causality which 
expresses itself in an explicit time ordering of current op- 
erators. It is very tempting to derive expressions which 
do not contain the time ordering any more. This can 
indeed be done but only for the first two cumulants [27l |. 

Thanks to intcgrability however, we can rely on yet 
another approach. Using the picture of single particle 
scattering, and extending the result for non interacting 
electrons |3 as has been done in the past for current and 
noise, we can identify xW as the average of product of 
single integrable quasiparticle events [(anti)kink bounc- 
ing off the boundary; we suppress the trivial prcfactor 



r/h from now on], 

x(A) = ( n {1 + [(e^' - - ^+) 

k 

+(e-^^-l)r7+(l-r7-)]}), (4) 

where the ri± are occupation numbers (equal to zero or 
one) for single particle kink and antikink states of mo- 
mentum k in the equilibrium (and unfolded) integrable 
model. The brackets denote average over the equilib- 
rium distribution determined by the external potential. 
Of course, this generating function correctly reproduces 
the known noninteracting results upon r]± = rjn/i^. We 
claim that, by a simple extension of the arguments used 
in [2^ I23 , it does give the correct result in the integrable 
case as well. Note that the expression (0J is rather im- 
plicit: there are in fact correlations between occupation 
numbers, which have to be correctly taken into account 
through the BA, and the product over k itself hides a 
complex integral, where the density of allowed levels de- 
pends on all the filling fractions. It seems possible to 
draw results for x(A) nevertheless, but here we will con- 
tend ourselves with using it to evaluate some cumulants. 
Defining N± = 77_|_(1 — ?7_) ± 77_(1 — 77+) and calculating 
the first three cumulants from C„ = (i)"9" lnx(A)/9A", 

C^ = Y^T{k){N^{k)), (5) 
fe 

C2=Y. T{k,)T{k2){N_{k,)N_{k2))c 

ki ,k2 

+ Y.T{k)[l-T{k)]{N+{k)), (6) 

k 

C3 = J2^ik) [1 - T{k)] [1 - 2T{k)] {N^{k)) (7) 
fe 

+ 3 T{ki)T{k2) [1 - T{k2)] {N-{ki)N+{k2))c 
+ T{ki)T{k2)T{k^){N^{ki)N^{k2)N^{k^)),, 

k^^k^^k^ 

where the subscript c denotes the normal ordering. Com- 
paring lO to lO and © to Eq. (27) of one imme- 
diately recovers the correct formulas for the current and 
the noise, while the third cumulant is given by the fol- 
lowing expression, which, along with Q), is the central 
result of the Letter, 



C3 — Cc + C/ + Cn , 



(8) 



where there is a correlation contribution just as in the 
case of noise, 



Cc = i^J J j dede'd9"r{e)T{6')T{e" 



X D{9, 



9', 9"), 



(9) 



with the correlation function D{9,9' ,9") = (A(P+ 
P_)(6i)A(P+ - P_)(6i')A(P+ - P-){9")) (and AP 



3 



P — P). The second ingredient is the current like contri- 
bution. 



cording to 



Ci^ j deT{6)[l-T{6)][l-2T{e)] 

x[P+(0)-P_(0)], (10) 

which indeed resembles the current Q up to the modified 
transmission coefficient. Finally there is a component 
which has the look of the second part in the noise, and 
which we shall call the noise like contribution, 

CN = 3eJ J d0d0'T{0)T{0')[l-T{d')] (11) 

x([P+(0)-P_(0)] {P+{0')[l-f_i0')] 
+ P^i0')[l-U{0')]}) . 

In all these formulas, £ is the length of the system (taken 
to infinity at the end), and brackets denote averages with 
respect to the equilibrium distribution in the presence of 
chemical potentials as determined by the TBA equations. 

The main problem now is to evaluate, using the TBA, 
the correlation functions in these expressions, and inte- 
grate over them. This is an a priori difficult endeavour, 
which is made possible by the use of several tricks gen- 
eralizing ideas used first in (2^. The easiest part is the 
current like contribution, which only requires the expec- 
tation values of level occupation probabilities P± : these 
come out from the numerical solution of the TBA equa- 
tions. The correlation component Cc can be extracted 
by an extension of the approach used in to evalu- 
ate D{0,0'). One first defines new chemical potentials, 
fi+{e) = -M-(^) = V/2T+xT{0) . Solving the BA equa- 
tions for quasienergies and calculating the free energy per 
unit length F one recovers Cc by a triple derivative ac- 



J 



Cc 



1 d^F 



T dx^ 



x=0 



(12) 



The background of this formula is the expression for the 
system's partition function as a functional integral over 
the level occupation probability, see the equation between 
Eqs. (5) and (6) of where the chemical potentials 
can be explicitly rapidity dependent. Finally, in order to 
calculate the integral in the noise like part Cat, one needs 
to introduce a new energy (per unit length) functional 
containing two c-numbers x and y as parameters. 



E = j de[[e+-ax{l- f^)+by]{0)P+{0) (13) 
- axil - /+) - bymP^{0) + ^ e,P,] , 



where the very last sum is over the breathers' degrees of 
freedom only and a{9) and b{6) are arbitrary functions of 
rapidity. Modifying the energy in this way and perform- 
ing the TBA program we find following equations for the 
quasienergies. 



T 



- 61') ln(l -he^'^+'^O 



- ^ jde' [$,+(0 - 0') + <i>,_(0 - e')] u{0')f-{0') , 

with e'^iCi = e± - a{0)x ± b{e)y + (a{0)xf^ for the 
(anti)kinks and e '- = e-,- for the breathers. The kernels <i>y 
can be found in |22 . l24j . The corresponding free energy 
is then 



F = 



-tj: d0 



e. In (1 + e-^'+^-) [e, - e[{l)/T + {1 + e^'"^') In (l + e-^-+''')] 



e;;(2)/r+ (1 + e^'-'^')ln(l + e-^'+A"-) + [a{e)x/T] J dO' [$,+ (0 - 0') + $,_(0 - 0')] f+{0')f-{9') 

I 



On the other hand, the same free energy can be found 
from the system's partition function |25l| . After subse- 
quent derivation with respect to both x and y, setting 
them to zero afterwards and fixing a{0) ~ T{0)[\~T{0)], 
b{0) ~ T{0) we immediately recognize the noise like con- 
tribution inj. 



C 



N 



3T 



d^F 



dxdy 



(14) 



One of the nontrivial checks of the formula (jSJ is the 
case V = 1/2, when there are no breathers . The results 
of [THlll^ are confirmed by a straightforward calculation. 
The second test is the accurate reproduction of the T = 



result, previously found in fl3 |. As can be immediately 
seen from the numerical calculation of the third cumulant 
as a function of temperature. Fig. ^ all correlation con- 
tributions vanish in the limit of low temperature. This is 
quite natural as in this case all thermal fluctuations die 
out. The only surviving contribution is then the current 
like part, which can be computed exactly in the same 
way as presented in leading, of course, to the same 
result. Many other checks are possible, and the structure 
of C3 is in fact quite fascinating - we will report about 
this elsewhere. 

Now, we would like to turn to the practical question of 
the measurability of the fractional charge through the 
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FIG. 1: Temperature dependence of the third cumulant and 
its constituents for V/Tb = 1. 
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FIG. 2: Temperature dependence of the square root of the 
ratio of the third cumulant and of current like contribution 
absolute values (IC3I and |C/|) to backscattered current Ibs 
as well as of the shot noise CI, to Ibs a.t v = 1/3 and Tb/V = 
0.1. The cusp is unphysical and comes about as a result of 
the Ci sign change. 



third cumulant. In the weak back-scattering limit at 
T/V = 0, it can be easily verified analytically that, 
when V/Tb ^ 1, the Schottky formula is valid C2/IBS = 
[here Ibs is the backscattered 



O 



{Tb/V) 



2(1-1.) 



current defined as Ibs = I{V,Tb = 0) - IiV,TB)], 
and similarly for the third cumulant C^/Ibs ~ + 

O {Tb /V)^''^~^^ . Looking at the detailed coefficients 



in these expansions shows that C2 gives slightly better 
results for the fractional charge at T = 0. At finite 
temperatures one has to use instead of the above C2 
only its shot noise constituent remaining after the sub- 
traction of the equilibrium current fluctuations: C2 = 
C2{V,T,Tb) - C2{V = 0,T,Ts). On the other hand, 
since the equilibrium fluctuations are Gaussian, the cor- 
responding third cumulant is identically zero, and no 
subtraction is necessary for the analysis of C3. That 
circumstance led Levitov and Reznikov to suggest that 
the measurement of the bare third cumulant would be a 
better way to access the value of v These authors 

also argued qualitatively that C3 should have a weaker 
dependence on temperature in the presence of generic in- 
teractions. We would like to point out though, that their 
discussion has been based on weak tunneling results for 
noninteracting particles. 

We are able to confirm this hypothesis precisely in the 
edge state tunneling problem. We find for instance that 
for a small value of Tb/V, while the shot noise gives a 
measure of the Laughlin quasiparticles charge only in the 
limit T/V small, the third cumulant has essentially no 
dependence on temperature over a huge range of values. 
This is illustrated in Fig. [21 A similar result holds if we 
keep T/V fixed and small, and vary Tb (i. e. the strength 
of the backscattering potential). Interestingly, such a 
remarkable behaviour of the third cumulant is entirely 
due to its current like component: in Fig. |21we plotted 
separately the fractional charge value induced by that 
contribution only. It appears to be quite insensitive to 
temperature change, only slightly oscillating just below 
the asymptotic v. 

To conclude, we have presented a procedure for ob- 
taining the cumulant generating function of the proba- 
bility distribution function to transfer a fixed amount of 
charge through a point contact in a FQH edge state de- 
vice. Using this approach we analyzed the behaviour of 
the third cumulant in detail. We predicted that the frac- 
tional charge can be much more reliably obtained from 
C3 than from the noise spectra and we identified the con- 
tribution to C3 which is responsible for this remarkable 
property. 

We would like to thank H. Grabert, A. O. Gogolin and 
B. Trauzettel for discussions. AK is Feodor Lynen fellow 
of the Alexander von Humboldt foundation (Germany). 



[1] L. S. Levitov and G. B. Lesovik, JETP Lett. 58, 230 
(1993). 

[2] L. S. Levitov, W. W. Lee, and G. B. Lesovik, Journ. 

Math. Phys. 37, 4845 (1996). 
[3] B. Reulet, J. Senzier, and D. E. Prober, Phys. Rev. Lett. 

91, 196601 (2003). 
[4] M. Kindermann and Yu. V. Nazarov, Phys. Rev. Lett. 



91, 136802 (2003). 
[5] W. Schottky, Ann. Phys. (Leipzig) 57, 541 (1918). 
[6] G. B. Lesovik, JETP Lett. 49, 592 (1989). 
[7] M. Biittiker, Phys. Rev. Lett. 65, 2901 (1990). 
[8] T. Martin and R. Landauer, Phys. Rev. B 45, 1742 

(1992). 

[9] R. de Picciotto, M. Reznikov, M. Heiblum, V. Umansky, 



5 



G. Bunin, and D. Mahalu, Nature (London) 389, 162 
(1997). 

[10] L. Saminadayar, D. C. Glattli, Y. Jin, and B. Etienne, 

Phys. Rev. Lett. 79, 2526 (1997). 
[11] L. S. Levitov and M. Reznikov, Phys. Rev. B 70, 115305 

(2004). 

[12] D. A. Bagrets and Yu. V. Nazarov, Phys. Rev. B 67, 

085316 (2003). 
[13] A. O. Gogolin and A. Komnik, Phys. Rev. B (2006). 
[14] H. Saleur and U. Weiss, Phys. Rev. B 63, 201302(R) 

(2001). 

[15] A. V. Andreev and E. G. Mishchenko, Phys. Rev. B 64, 
233316 (2001). 

[16] M. Kindermann and B. Trauzettel, Phys. Rev. Lett. 94, 
166803 (2005). 

[17] A. Komnik and A. O. Gogohn, Phys. Rev. Lett. 94, 
216601 (2005). 

[18] L Safi and H. Saleur, Phys. Rev. Lett. 93, 126602 (2004). 
[19] X. G. Wen, Phys. Rev. B 41, 12838 (1990). 



[20] X. G. Wen, Phys. Rev. B 43, 11025 (1991). 
[21] C. L. Kane and M. P. A. Fisher, Phys. Rev. B 46, 15233 
(1992). 

[22] P. Fendley, A. W. W. Ludwig, and H. Saleur, Phys. Rev. 

B 52, 8934 (1995). 
[23] S. Ghoshal and A. B. Zamolodchikov, Int. J. Mod. Phys. 

A 9, 3841 (1994). 
[24] P. Fendley, A. W. W. Ludwig, and H. Saleur, Phys. Rev. 

Lett. 74, 3005 (1995). 
[25] P. Fendley and H. Saleur, Phys. Rev. B 54, 10845 (1996). 
[26] Yu. V. Nazarov, Ann. Phys. (Leipzig) 8 (SI-193), 507 

(1999). 

[27] G. B. Lesovik and N. M. Chtchelkachev, JETP Letters 
77, 393 (2003). 

[28] In all cumulants we omit the trivial prefactor e" as in ac- 
cordance with our definition of P{Q) we are counting the 
number of elementary charges and not the actual charge 
amount. 



